Interplay of disorder and nonlinearity in Klein-Gordon models: Immobile kinks 
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We consider Klein-Gordon models with a 5-correlated spatial disorder. We show that the prop- 
erties of immobile kinks exhibit strong dependence on the assumptions as to their statistical dis- 
tribution over the minima of the effective random potential. Namely, there exists a crossover from 
monotonically increasing (when a kink occupies the deepest potential well) to the non-monotonic 
(at equiprobable distribution of kinks over the potential minima) dependence of the average kink 
width as a function of the disorder intensity. We show also that the same crossover may take place 
with changing size of the system. 
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I. INTRODUCTION 

An extensive research work on the static and trans- 
port properties of nonlinear excitations in various soliton- 
bearing disordered systems has been undertaken in the 
last decade (see Refs. and 0). It is well known that 
when taken separately both the nonlinearity and disor- 
der contribute to the localization effects, the character of 
localization being however essentially different. The non- 
linearity results in the possibility of existence of nonlinear 
localized excitations usually referred to as solitons that 
are rather robust and can propagate through the system 
undistortedly. At the same time the disorder (in linear 
systems) evokes the Anderson localization, which mani- 
fests itself in the behavior of the transmission coefficient 
of a plane wave decaying exponentially with the system 
width. These two localization mechanisms are competi- 
tive to some extent; taken together, the nonlinearity and 
disorder may lead to a number of qualitatively new ef- 
fects, namely, the transmission coefficient tending to zero 
with.Jncreasing the system length according to a power 
lawtfQ rather than exponentially (which would be oth- 
erwise the prop|erJtjy of a linear system); there can arise 
a multistabilityu'cl in the wave transmission through a 
disordered slab; excitations in highly nonlinear or multi- 
dimensional nonlinear Schrodinger (NLS) systems (which 
would either dispe]fse-|Or collapse otherwise) can be sta- 
bilized by disorder .I3i3 

In the present paper we study the static properties of a 
one-kink solution (or, equivalently, of a diluted kink gas) 
of disordered Klein-Gordon models where the disorder 
is assumed to be a ^-correlated Gaussian spatial noise. 
The disorder of the kind is akin, for example, to Joseph- 
son junctions, where it is caused by the fluctuations of 
the gap between two superconductor plates. Quite re- 



cently Mints proposedEi a similar model with randomly 
alternating critical current densitji-to account for a self- 
generated magnetic flux observedli^ experimentally. We 
discuss his results in more detail in the Conclusion. 

In gengj^, Klein-Gordon models have been repeatedly 
studiedElalld for the disorder represented by a lattice of 
(5-like impurity potentials with random positions of the 
impurities and either equal or randomly distributed in- 
tensities. It was found out that in a number of cases 
the kink dynamics in disordered systems could be ad- 
equately described within the framework of the collec- 
tive coordinates approach (Refs. [l4-[l5|, and references 
therein). This is the background of our restricting our- 
selves to the said approach in the scope of the present 
paper providing however at each step the validation of 
analytical results comparing them to the numerical sim- 
ulations of the original system. We emphasize that in 
opposite to, for example, Refs.. B and [ll| , we use Rice's 
collective coordinates approachOO with the kink width 
being the variational parameter. 

A similar approach has been recently applied for 
two other cme-dimensional (ID) systems: Bussac et al. 
investigatedB the effects of the polaron ground state in g 
deformable chain, while Christiansen et al. consideredQ 
the stabilization of nonlinear excitations by disorder in 
the NLS model. It should be indicated that investigation 
of these two (in fact closely related) systems leads to 
one and the same result: the width of stationary solitons 
decreases with growing intensity of the disorder. The im- 
portance of this conclusion resides in its prediction that 
the disorder can stabilize otherwise unstable solitons in 
2D and 3D NLS models. Qjiite recently this prediction 
was borne out numericallya for the 2D case. It must 
be emphasized that for NLS models the conclusion does 
not depend on the averaging procedure: one can equally 



perform averaging either on absolute ground statesQ or 
over all local mininiA, of the effective random potential 
with equal weights J3'Q We show that it is not the case 
for the Klein-Gordon models: their properties exhibit 
strong dependence on the assumptions as to the statis- 
tical distribution of kinks over the minima of the effec- 
tive random potential. For the purely dynamical problem 
these statistics are left beyond consideration and should 
be thus imposed as an additional assumption. We use 
Jaynes's maximum entropy inferencalj for this purpose. 
The outline of the paper is the following. In Sec. 
we present the model and derive the equations for the 
collective coordinates of the kink taking into account the 



III 



disorder via the effective random potential. In Sec 
we investigate both analytically and numerically the case 
of immobile kinks and demonstrate the existence of the 
crossover between monotonic and nonmonotonic depen- 
dence of the average kink width as a function of the dis- 
order intensity. In Sec. fV we summarize the exposed 
results. 



II. COLLECTIVE COORDINATES APPROACH 

We consider a Klein-Gordon (KG) model in the pres- 
ence of space disorder. The Hamiltonian of the system 
has the form 

H = I" dxl^^if, + (fl) + [1 - ^e(x)]<I>(^)| , (1) 

where the subscripts stand for partial derivatives with 
respect to the indicated variables and units are chosen so 
that the Hamiltonian is already in the scaled form. The 
potential $(</)) has the form 



$((/)) = 1 — COS0 

for the sine-Gordon (SG) model, and 



(2) 



(3) 



for the (f)^ model. We assume that e(a;) is delta-correlated 
spatial disorder 



{e{x)e{x')) = 6{x - x') , 



(4) 



(the brackets (...) denote averaging over all realizations 
of the disorder) with the Gaussian distribution 



effects of disorder in both these models. But since the 
qualitative features of the exposed models turned out to 
coincide in the scope of our present investigations, we 
dare not overload the paper with unnecessary repetitions 
and restrict ourselves with presenting in detail the sine- 
Gordon model only, keeping in mind although that every 
stage of the calculations applies for the 0* model as well. 
The SG system is governed by the equation of motion 



(/>** - 0XX + [1 - V^ix)] sin + 7(/)t = , 



(6) 



where the damping term with the damping constant 7 
has been included. It is well known that in the absence of 
disorder and damping ('7 = 7 = 0) Eq. (o) is completely 
integrable and possesses a topologically stable solution in 
the form of a kink given by 



(pK {x,t) —A arctan exp 



x-X{t) 

m 



(7) 



where X{t) = Xq + vt is the kink coordinate, v is its 
velocity, and L = \/l — v^ is the kink width. 

In the general case of Eq. (13) for a numbp; of situations 
the kink emission is exponentially smalljid so that the 
kink dynamics can be studied by the collective coordinate 
approach. In the framework of this approach the vari- 
ables X(t) and i(i) are understood as time-dependent 
variational parameters. Inserting Eq. (0) into Hamilto- 
nian (|l|) as a trial function, we obtain the effective Hamil- 
tonian 

^e// = ^v\ + ^L + V{L) + y({e}, L, X) , (8) 
with momenta 



8 dX 



_ 2TT^dL 



Here 



C/(L) = -+4L 

1j 



(9) 



(10) 



is the potential function in the case of no disorder and 



/oo 
dxe{x) sech^ 
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x-X 
L 



(11) 



is the effective random potential arising because of the 
disorder term. Then, taking into account the .damping, 
we arrive at the following equations of motionita 



dp_ 



X 



dt 



IPX 



dX 



y({e},L,X) = 



(12) 



pH^)] = ^cxp[-e2(a;) 



We have studied both SG (|) and (j}^ (|) models. Al- 
though properties of these two models show many simi- 
larities, they also exhibit a number of interesting distinc- 
tions related, in particular, to the existence of a breathfii; 
state in the SG model and of Rice's internal modelij'tZl 
in the (j)^ model. So, it was important to compare the 



dpL 
dt 



IPL 



(5) ^+.„.+M + ^ 
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^[^(^) 



F({e},L,X)] = 0. 



(13) 



In the following section we solve approximately these 
equations of motion for immobile kinks and compare the 
results to the results of direct numerical integration of 
Eq. d). 



III. RESULTS FOR IMMOBILE KINKS 



V{{e},L,X)^27^^,{{e},X) 



(21) 



As a consequence of the damping 7 the kink wiU even- 
tuaUy stop at some stable or metastable stationary posi- 
tion along the system. Here we do not consider this tran- 
sient stage and assume that the kink is already immobile 
[dpL/dt = PL = Px =0]. In this case the equations of 
motion (n2h and (jl3|) take on the form 



d 
dX 



y({e},L,X) = 0, 



— [UiL) + Vi{e},L,X)] 



(14) 



(15) 



Considering the center-of-mass motion described by Eq. 
(O) we observe that for each realization of the ran- 
dom potential e{x) the stable stationary position X ~ 
Xm{{£},L) of the kink is defined by the point where 
V{{e}, L,X) has a minimum with respect to X. Thus 
we can now insert the value X = Xm({e},_L) into Eq. 
( p^ ) and, solving the resulting equation 



l + i^^«^>'^'^" 



-1/2 



(16) 



get the value of the stationary kink width L{{e}) at given 
position of the kink. Encountering a similar problem for 
the case of NLS system, Christiansen et al. invokedQ the 
mean-field approximation 

{A.Vi{e},L,X,^)) « ±{Vi{e},L,X^)) . (17) 

Further the estimation of the quantity|_/t^) was per- 
formed using Rice's averaging theoroBLEj-EH However, be- 
ing rather good for the NLS modclQ the mean-field ap- 
proximation (^7|) fails for the KG models. Thereby we 
were forced to use a more precise averaging procedure 
calculating (dV/dL) directly. 

Expanding Eq. (H^) into series up to the second order 
in 77, 



IdV 3 fdV 



(18) 



and solving by iterations we get after averaging [which 
by means of Eq. (0) can be performed for the terms 
containing 77^ exactly] the average kink width 



L,ar^{L{{e}))^l + j{X{{e},X„,)) 
where the averaging of the function 
\{{e},X)- 



180 



dxe(.)(x-X)i^^^^p^ 
cosh (a; — X) 



(19) 



(20) 



is performed over all realizations of the disorder in the 
points X = X„i in which the potential 



with 



/C30 
dxe{x)sech^{x-X) (22) 

-00 



takes on its minima on X . 

Thus we arrive at the problem of performing the av- 
erage of A({e}, Xm) over the minima Xm of the function 
/i({e},X). It is convenient for later use to perform this 
averaging in two steps calculating at the outset 



A{p) 



XPiiX\p)dX 



(23) 



and thereafter 



{X{{e},X^))= / A{fi)P„,{fi)dfi . 



(24) 



Here P/(A | jl) is the conditional probability that 
X{{e},Xm) has the value A if p{{e},Xm) equals to /i. 
Correspondingly, A(/2) is the value of X{{e},Xm) av- 
eraged over all realizations of the disorder for which 
^{{e},X„i) is equal to p. It is difficult to calculate A(u) 
analytically but numerical simulations show (see Fig. m 
that up to very good accuracy the dependence A(/i) is 
linear: 



A(/i)~ -0.344/2-0.500 



(25) 
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FIG. 1. The dependence A(/i) found numerically. The cal- 
culations are performed for a million minima of the function 
pi{e},X). 



Substituting it into Eq. (|2J) we obtain that 
(A({e}, Xra)) ^ -0.344 {^l{{€}, X,^)) - 0.500 , (26) 
where 



(Ai({e},X„0) 



ilPra{iJ)dfl 



(27) 



is the average value of the function ii{{e},X) over its 
minima X^- Here the probabihty density Pmifi) is a 
product of two factors. The first one is the probabihty 
density that some arbitrary chosen minimum Xm of the 
function /i({e},X) will be equal to jl. We denote this 
probability density as PminiiJ-)- The second factor is the 
conditional probability that if the minimum is equal to 
jl it will be actually occupied by the kink. It is evident 
that in a real system the kink is more likely to occupy 
the deeper minimum than the shallow one. So to be 
consequential one must ascribe to every minimum of the 
function iJ,{{e},X) some probability weight and average 
taking into account those probabilities. But the values 
of these probabilities are in general determined by the 
whole prehistory of the kink. These values are not con- 
tained in the dynamical equations of motion that state 
only that the kink should take on some minimum regard- 
less to its depth. It would be a cumbersome problem to 
calculate them appropriately. That is why two limit cases 
are in general ccmsidered: either the kink seats itself into 
the deepest wellQ ot. it rather occupies any of them with 
equal probabilities. Q As it was already remarked in the 
Introduction, both limit cases lead to qualitatively the 
same results for the NLS model. 

It is not the case of KG models where, as it will be 
shown later, different assumptions as to the a priori 
weights of minima lead to qualitatively different behavior 
of the kink width. Since we are merely lacking informa- 
tion sufficient enough to reconstruct these weights in an 
objectivistic fashion, axemedy would be Jaynes's maxi- 
mum entropy inference, tj according to which the simplest 
self-consistent unbiased choice is to assume that the kink 
will occupy a potential well corresponding to the mini- 
mum Xm of the function p,{{e}, X) [for given profile e{x)] 
with probability proportional to e"^'^'^^'^^''^'"^ thus intro- 
ducing an additional parameter (3 (following Jaynes we 
shall call it a conjugate parameter). By this expedience 
we present some natural interpolation covering two men- 
tioned limit cases: of equiprobable distribution (/3 — 0) 
and of averaging over the deepest minima only (/3 -^ oo). 
So we can write 



1 



-Pt^r 



where 



-^m vMJ ^^ Pmin 



if^), 



Z{f3)= / rfAe-'^>™„(A) 



(28) 



(29) 



To calculate Pminip-) we fcjUoss^ Ref. and make use of 
the Rice's averaging theorerroGj [valid for the case, well 
attested by the numerics, of /i({e}, X) being a stationary 
centered Gaussian process] stating that the probability 
density of some given minimum of the function ^({e}, X) 
to be equal to p, is 



(30) 



1 ( p- 


J I ^"- 


^"""^^^ V2^FMo"lV2]Wo' 


\ MoMi^ 


where the function 




cr(y, k) = Ke"(^/'')' 




- 2v/l - «2ye-^^ / 


e-'"dt 



l(y/K)Vr- 



(31) 



and the spectral momenta 



Mo = {H{e},X)]^) = - 

Ifi 

M,^{[^^x{{e},X)f) = - 

64 
21 



(32) 



M^^([^lxx{{e},X)]') 



were introduced. 

Hence the partition function takes on the form 



1 r 



dxe 



8/3/32; 



a I x,-\^ 



(33) 



and can be expanded into series in /3 yielding the average 
value 



{^,({e},X,n))=-—\nZiP) 
'y-^(64-7^)/3-KO(/32). 



(34) 



And now, substituting it into Eqs. 
obtain for the average kink width 

Lvar- 1 + (0.193/3 -0.059) 77 



and (|l|) 



we 



180 



V 



(35) 



plays part of the partition function. 



Thus it is seen that there is a qualitative change of the 
kink width behavior as function of the disorder intensity 
77 according to whether the value of conjugate parameter 
/3 is below or above some critical value Per ~ 0.3. 

At small f3 the average kink width is a nonmonotonic 
function of the disorder intensity rj: it decreases at small 
intensities but starts to increase thereafter. This result is 
well attested by the direct numerical calculations of sta- 
tionary kink solutions of the initial equation of motion 
(0). In Fig. we compare the analytical prediction given 
by Eq. ( pq ) with the numerical results for the case of 
equiprobable distribution of the kinks over the potential 
wells {P — 0). The numerical results have been obtained 



as an average of 1000 realizations of the disorder. Two 
different expressions for the kink width were calculated: 



Lcos == 7 / ^2;{1 - cos[(l){x)]} , 



and 



Lder — 8 



dx 



d(f){x) 



dx 



(36) 



(37) 



From the point of view of the collective coordinate ap- 
proach Lcos and Ld^r should coincide with the value of 
L introduced in Eq. (0). Indeed, as it is seen from the 
figures, it does take place for < r/ <^ 0.2; thus, these 
are limits where the collective coordinate approach works 
well. 
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FIG. 2. The dependence of the average (on 1000 realiza- 
tions of the disorder) kink width {L{{e})) vs disorder intensity 
77 under the assumption that the kink occupies every potential 
well created by disorder with equal probability. The system 
length i? = 15. 




FIG. 3. The dependence of the average (on 1000 realiza- 
tions of the disorder) kink width (L({e})) vs disorder intensity 
77 under the assumption that the kink sits into the deepest po- 
tential well created by disorder. The system length R — 15. 

Going on with the analysis of Eq. (^) we see that at 
big values of the conjugate parameter (/3 > (3cr) when 
the kink rather occupies deep potential wells, the kink 



width should grow monotonically with the disorder. But 
it is evident that in this case the analytical approach dis- 
cussed above is applicable to systems of infinite length 
only. For the finite system the case of big /3 represents 
the situation when the kink sits in the deepest potential 
well. Obviously its average depth essentially depends on 
the system length. We can make an. estimation of this de- 
pendence drawing on the formulacj for the average num- 
ber of minima of the function /i({e}, X) on an interval R 
whose values lie below some fl: 






(38) 



Inserting there Nmm = 2 (/i is an absolute minimum on 
the interval R but not on the longer interval) one can 
estimate its average value 



{K{^}^Xabs.rn)) 



and, substituting it into Eqs. (uM and ( p6D one can find 
that the average kink width equals 




"■^'■""Isfcl-^-^^l" 
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(40) 
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FIG. 4. The dependence of the average (on 1000 realiza- 
tions of the disorder) kink width {L{{e})) vs disorder intensity 
■q under the assumption that the kink sits into the deepest po- 
tential well created by disorder. The system length R = 5. 

It is seen that for finite-size systems the character of 
dependence L^ariv) depends on the size of the system 
R. It is interesting that even for the case of averaging 
over the absolute minima considered here, the function 
Lvarivi) grows monotonically with rj only for the systems 
that are large enough {R > 7.5). The reason is that for a 
small system the number of potential wells of the effective 
random potential is too small to yield the average over 
absolute minimum that would be essentially smaller than 



the average value calculated over all minima. Indeed, 
Figs. H and ||, in which we compare Eq. ( [JOJ ) to the re- 
sults of the numerical calculations for i? = 15 and R = 5, 
lend support to the validity of the approach leading to 
Eq. (M) . One can see from these figures that the average 
kink width grows monotonically with rj for R = 15 but 
is nonmonotonic (similar to the case depicted on Fig. ||) 
for small R — 5. But in this latter case the boundary 
conditions become very important and most likely the 
are responsible for the difference between Figs. 



IV. CONCLUSION 



cely tney 
I and I 



al He has shown that in the case of 77 < 1 a state with a 
self-generated flux exists and can be studied experimen- 
tally in the presence of Josephson vortices. However, as is 
shown in the present paper, the Josephson energy [equal 
to 4:Lcos in Eq. (pq)] and the magnetic energy [equal to 
4/ Lder in Eq. (|3J)] of the Josephson vortices are func- 
tions of the intensity of fluctuations of the insulating layer 
thickness. And their contribution into the experimen- 
tally observable magnetic flux will strongly (up to the 
change of a sign) depend on the statistical distribution 
of the vortices along the Josephson junction. Thus, for 
the proper description of the problem one must develop 
a thermodynamic model. 



In the paper we consider the Klein-Gordon models with 
the ^-correlated spatial disorder and investigate both an- 
alytically and numerically the width of immobile kinks 
as a function of the intensity of disorder. The analytical 
collective coordinates approach is based on Rice's ^'•Wn 
aging theorem from the theory of random processestj'E2l 
as well as .oii the maximum entropy inference proposed 
by Jaynes.tSl 

We have shown that the properties of the kinks exhibit 
strong dependence on the assumptions as to their statis- 
tical distribution over the minima of the effective random 
potential. Namely, there exists a crossover from mono- 
tonically increasing (when a kink occupies the deepest 
potential well) to the nonmonotonic (at equiprobable dis- 
tribution of kinks over the potential minima) dependence 
of the average kink width as a function of the disorder in- 
tensity. We have shown also that the same crossover may 
take place with the changing size of the system: the av- 
erage kink width monotonically increases for the systems 
of big size but is nonmonotonic for the small ones. 

It is interesting to compare the effects of the disor- 
der in the KG model with the effects in the nonlinear 
Schrodinger (NLS) model. As it was recently shown in 
Refs. B-H, the 5-correlated spatial disorder in the NLS 
systems creates an additional factor contributing to the 
decrease of the excitation width. This effect, being in- 
sensitive to the manner of the statistical distribution of 
kinks over the minima of effective random potential, fa- 
vors the stabilization of excitations in highly nonlinear or 
multidimensional systems, which would either disperse or 
collapse otherwise. The stabilizing function of disorder 
is of no doubt important for practical applications and 
to elucidate the extent to which it is universal seems to 
be an intriguing question. The considered example of 
the KG systems demonstrates that there exists a class of 
systems, for which, in contrast to the NLS system, the 
effects of disorder can lead in different cases to diametri- 
cally opposed behavior. 

In the case of the SG model we can consider the term 
r]e{x) as a change of the Josephson current density due to 
fluctuations of the thickness of an insulating layer. Quite 
recently Mints studiedEI such a model-|to account for a 
self-generated magnetic flux observedll^ by Mannhart et 
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